Quantum disorder versus order-out-of-disorder in the Kugel-Khomskii model 
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The Kugel-Khomskii model, the simplest model for orbital degenerate magnetic insulators, ex- 
hibits a zero temperature degeneracy in the classical limit which could cause genuine quantum 
disorder. Khaliullin and Oudovenko [Phys. Rev. B 56, R14 243 (1997)] suggested recently that 
instead a particular classical state could be stabilized by quantum fluctuations. Here we compare 
their approach with standard random phase approximation and show that it strongly underestimates 
the strength of the quantum fluctuations, shedding doubts on the survival of any classical state. 



Motivated by the developments in the manganites, the 
interest in the role of orbital degeneracy in strongly cor- 
related systems revived. A classic model in this context 
is the Kugel-Khomskii model , believed to be realistic 
for KCuFs and related systems (one hole per site, de- 
generacy of the Cg orbitals). We recently discovered that 
this model poses a rather fundamental problem ||^ : in the 
classical limit a point exists in the space of physical pa- 
rameters where its ground state becomes infinitely degen- 
erate, due to a novel dynamical frustration mechanism. 
This classical degeneracy is lifted on the quantum level, 
and by analyzing valence-bond type variational states 
we arrived at the suggestion that the ground state for 
5 = 1/2 might well be an incompressible spin fluid. In a 
follow up paper, Khaliullin and Oudovenko [Q suggested 
that instead the quantum fluctuations act to single out 
one particular classical state (the one with Noel order and 
d^y2_j.2 orbitals occupied by holes) over all others by an 
order-out-of disorder mechanism: the classical degener- 
acy is lifted by the differing strength of the fluctuations 
around the various classical states but these fluctuations 
are not severe enough to destroy the classical Neel or- 
der completely. Their suggestion was based on a par- 
ticular decoupling scheme and in this communication we 
will demonstrate that for rather simple reasons this de- 
coupling scheme implies a serious underestimation of the 
strength of the fluctuations, shedding serious doubts on 
the possibility that classical order survives after all. 

The Kugel-Khomskii model describes a three-dimen- 
sional (3D) cubic Mott-Hubbard insulator with a single 
hole/electron in Cg orbitals {x^ — y^ '^^^ l^;), 3z^ — l \z)), 
possessing, in the absence of virtual hoppings, orbital 
degeneracy in addition to the standard spin degeneracy. 
Its minimal version is given by {J — jU being the an- 
tiferromagnetic (AF) superexchange, with t the hopping 
element between \z) orbitals along the c-axis) 



if = J^[4(5..5,)(rf-i)(r;-i) 



(1) 



where Ez is the energy splitting between the Cg orbitals, 
acting as a "magnetic field" for the orbital pseudo-spins. 
It is used to investigate the system when it approaches 
the degeneracy point — 0. The spin operators Si are 
S — 1/2 spins, while the orbital degrees of freedom are 
described by (2 x 2) matrices in the pseudospin space, 



a{b) 



(2) 



(?)■ 



and a selects the cubic axis (a, 6, or c) that corresponds to 
the orientation of the bond (ij). The cr's are Pauli matri- 
ces acting on the orbital pseudo-spins ja^) = o ) ' 1-^) ^ 
^ 1 . Hence, a Heisenberg model for the spins is coupled 

into an orbital problem. Here we ignore the (physically 
important) niultiplet splittings due to a finite value of 
the atomic Hund's rule coupling (Jh), and focus on the 
special point Ez, Jh 0, contained in model Eq. (^: it 
is easy to see Q that in the classical limit the system is 
dynamically frustrated and an infinity of classical phases 
become degenerate at zero temperature. This degeneracy 
is lifted on the quantum level and one expects quantum 
effects to take over at this point itself, as well as in its di- 
rect vicinity ||], in analogy to what seems established in 
geometrically frustrated spin models |^ . If a disordered 
state would be stabilized by quantum effects, orbital de- 
generacy could be added to the list of mechanisms leading 
to a spin-liquid, such as the frustrated Ji — J2 Heisen- 
berg antiferromagnet (HAF) the bilayer HAF and 
two-dimensional (2D) lattices with a reduced number of 
magnetic bonds, as reahzed in CaV409 |^]. 
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Quite generally, the transverse modes may be 

calculated starting from the equations of motion |^], 

mS+\-)) - ...]) + (([S+,H]\...)), (3) 

- ^([^^ ■■•]> + mt,H]\...)), (4) 

and using a generalization of the LSW theory. Here 
is either or = S^cr^, while /C^ is cither 
or K+- = Sfar. The first pair of 



++ — 



S, 



Green functions stands for spin-wave (SW) excitations, 
while the second pair describes mixed spin- and- orbital- 
wave (SOW) excitations. Similarly a longitudinal mode 
is given by 
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E{{at\...)) ^ -{[af ,...]) + {{[at ,H]\...)), 



27r 



(5) 



where the Green function describes a purely orbital ex- 
citation. At each site the full set of local operators de- 
scribing these excitations constitutes a so(4) Lie algebra. 
Even though the spin- wave operators form a subalgebra, 
as seen from the familiar su{2) commutators together 
with the additional commutators 

\St,Stat]^-St, \St,Stat] = -St. (6) 

while the same holds for the spin-and-orbital operators, 

[Kt^.St] = -Kt^, [Kt^.K^^] = 45f ± 2a,f , (7) 

one needs with Hamiltonian (|^) for the calculation of the 
SW and SOW excitations also commutators like 



-25'^ 



(8) 



Glearly, the SOWs cannot be separated from the SWs, 
and one has to solve simultaneously Eqs. (||) and (^). 

The random-phase approximation (RPA) for spinlike 
operators linearizes the equations of motion by the famil- 
iar decoupling procedure 



{{AB,\...))^{A){{B,\...)) + {B,){{A.\...)). 



(9) 



It is crucial that the decoupled operators Ai and Bj have 
different site indices, so that this procedure docs not vio- 
late the local Lie-algebraic structure of the commutation 
rules (j^H). In the Neel-type AF phases with either \x) 
(AFxx) or \z) (AFzz) orbitals occupied, one now finds 
after Fourier transformation and using the nonzero ex- 
pectation values of 5*1, cr|, and 5*1 (t| operators, two ex- 
citations (a 

,(n)l2 



^ k 



X, z for AFxx and AFzz, respectively), 
i^l - rlf - 2{Xl 



rl)iQl,-Rl) 
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2P2 )2 



1/2 



(10) 



X — 1 

^x(z) 2 ' 



'x(z) — 2 -^^ 



The orbital dependence enters the fc-independent field, 

(11) 

with £2 = Ez/ J , and the dispersion is given by, 

Q.k = h+ik). = h+ik) + 47.(fc), (12) 

^xfe-|^^7-(fc). n£ = 5V37~(fc), (13) 
Rk = |7+(fc), (14) 

with 7±(fc) = 5 (cos kx ± cos ky) and ^z{k) — cos fc^. 

The dispersions of SW and SOW are shown in Fig. 
|l|. It is straightforward to verify that the SW disper- 
sion is 9J/2, given for the AFxx phase by the superex- 
change of 9 J/4 between \x) orbitals in the (a, 6)-planes, 
and for the AFzz phase by strong interactions of 4 J along 
the c-axis and weak superexchange of J/4 in the (a, b)- 
planes. In both phases one finds that the coupling be- 
tween the modes due to the P^j: ^ "f-{k) term is strong, 
and the excitations have pure character only in the planes 
of 7_(fc) = 0, as seen along T — L{K) lines. In particular, 
this coupling increases along the T — X direction, and 
precisely compensates the dispersion due to the orbital 



dynamics 7+(fc). This results in a soft mode w^^^ = 
along the F — X{Y) direction in both AF phases. As 
we have shown before j^, finite masses are found in the 
directions perpendicular to the soft mode lines, which 
gives a logarithmic divergence of the quantum correction 
to the order parameter, {5S^) ~ InA^, with Aj — > for 

KhaliuUin and Oudovenko Q calculate instead a SW 
and a longitudinal excitation {{at\...)) ^) first, and then 
include the effect of orbital fluctuations in the transverse 
channel (our SOW) in a perturbative way. Their selfcon- 
sistent calculation gives a finite energy and weak disper- 
sion of the orbital mode (in the present RPA approach 
the orbital excitation is found at w = 0). This approach 
violates the commutation relations in the Lie algebra (^ 
^) , and only for this reason the SW and SOW excitations 
become independent from each other. In the present RPA 
language it implies that composite spin-and-orbital oper- 
ators, Sfa^, are factorized into independent products of 
spin {Sf) and orbital (trf ) operators separately, and the 
commutators given by Eqs. (^-|^) effectively either van- 
ish, e.g. [St,Slaf] ^ Sf]{af) = 0, or give a differ- 
ent result, e.g. [St^Sfaf] ^ [5+,5f](af) = 
We call this procedure the SW-I-SOW scheme; it is for- 
mally equivalent to assuming P^j: = in Eq. (p^). The 
SW modes depend now solely on the actual magnetic in- 
teractions, while the SOW modes are identical in the two 
phases and the soft mode behaviour is absent (Fig. |^). 
However, the SOW mode is gapless at the F point as a 
consequence of the rotational invariance of the classical 
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AF order when the occupied orbitals are rotated between 
|a;) and \z). 

We calculated the order parameter (S"^) in both AF 
phases including quantum corrections using a generalized 
RPA approach which leads to the identity, 

(SDrpa {SrS+) - \ (SrarS+at). (15) 

where i £ A, and A is the "f-spin sublattice. The iden- 
tity (|l^) follows from the expansion of the Sf opera- 
tor in the so(4) algebra and replaces the su{2) relation 
(^Sf) = ^ — {S~S^), familiar from the Heisenberg model. 
It includes the renormalization due to both transverse 
modes in the spin-orbital model (|l|). Similarly the or- 
bital occupancy is renormalized by {cr^crf) fluctuations 
due to the longitudinal mode. The correlation functions 
are found from the respective Green functions 

= j du: I 1^23((%|^,,))^_„ j . (16) 

-oo \ fc / 

Eq. (|l|) reproduces the result for the 2D HAF, (Sf ) ~ 
0.303, in the limit of +oo, while {Sf) ~ 0.251 for 

the strongly anisotropic 3D HAF at Ez — > — oo. The 
values of {Sf) are, however, strongly reduced when the 
degeneracy point [E^ = 0) is approached (Fig. ||), and 
the quantum corrections overshoot the mean-field value 
(S'f )mf for -0.04 < Ez/J < 0.30, and diverge at E^ = 0. 
In contrast, these corrections are much reduced within 
the SW+SOW scheme, and the divergence at E^ is re- 
moved ((S'f) ~ 0.05 in both phases). This is qualitatively 
equivalent to the results of Ref. Q , where the renormal- 
ization of (5f ) due to the SOW was included only pertur- 
batively, and a value 0.191 was found in the AFzz phase. 
This somewhat smaller quantum correction results from 
the finite gap in the orbital excitation. 

Further evidence that the stability of the LRO phases 
is overstimated in Ref. Q comes from energy calcula- 
tions. For convenience we define the ground state energy 
per site as a quantum correction beyond the mean-field 
value, 

E=^{H)+Ez{tI) +?,,]. (17) 

A simple estimation aX E^ — Q using the Bethe ansatz so- 
lution for a disordered one-dimensional (ID) chain along 
the c-axis, and no magnetic correlations in the (a, b)- 
planes gives E = —0.648 J while a somewhat better 
energy of —0.656 J was obtained using plaquette valence 
bond (PVB) states either with singlets alternating along 
the a- and fe-axis in the (a,5)-planes (PVBA phase), or 
with single planes of such alternating singlets interlay- 
ered with two planes of singlets along the c-axis (PVBI 
phase), as explained in Ref. pj. 



For the LRO phases, in spite of the divergent correc- 
tion to the order parameter, an energy can still be es- 
timated using the RPA corrections for the symmetry- 
broken state. Here these estimates starting from the 
states with LRO give lower energies than the above sim- 
ple estimates for the VB states. This is not surprising, 
as it is known that improved VB wave functions that in- 
clude the resonance between spin singlets lead typically 
to large energy gains, but are difficult to treat already in 
spin models ||l^. First, within the generalized RPA ap- 
proach one finds the largest quantum corrections in the 
AFxx phase (Fig. This shows that the AFxx phase is 
'more unstable' against disorder, in agreement with intu- 
ition and with Ref. jj] . We believe that the lowest energy 
— 0.896J obtained in the AFxx phase at E'^ = comes 
close to the true ground state. This is consistent with the 
experience with the ID HAF, where one finds an energy 
of -0.429J using the LSW theory, which is only 3.2% 
above the exact value — 0.443J. We note that the en- 
ergy obtained within the simplified SW-I-SOW approach 
is much higher, even above that of the disordered phases 
(PVB states). In contrast, the SW-I-SOW approach gives 
for the AFzz phase an energy somewhat lower than that 
of the PVB states, and our value of E differs only by 
0.005 J from that reported by KhaliuUin and Oudovenko 
in their scheme [Table I]. This indicates the qualitative 
similarity of these two approximations in treating the 
quantum fluctuations related to simultaneous spin and 
orbital flips (SOW excitations) ~ in both cases the effect 
of such fluctuations is severely underestimated. 

Summarizing, the results presented in Ref. [Q are in- 
conclusive, as their approximation violates the so(4) dy- 
namical algebra describing the microscopic excitations. 
In contrast to the result of the perturbative treatment 
of Ref. 1^ , the RPA calculation yields an unstable AFzz 
(and AFxx) phase at orbital degeneracy, as also found in 
spin systems Q . As the LSW theory performs quite well 
in simple spin systems with S — 1/2 | [Tl[ |, this strongly 
suggests that the ground state of the Kugel-Khomskii 
model is a spin-liquid. To our knowledge, the present 
case is unique in the sense that singlets arranged in a 
3D valence bond solid (PVB states) allow a lower energy 
than that of a classical state. However, it might well 
be that the final verdict on these matters has to wait for 
the systematic approach to the quantization of classically 
frustrated problems, which is still to be invented. 
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TABLE I. Ground state energy E, in units of J, as ob- 
tained for the Kugel-Khomskii model using the full RPA 
(RPA) and decoupled SW and SOW excitations (SW-^SOW), 
compared with the energy found in Ref. Mj. 



method 


AFzz phase 


AFxx phase 


RPA 


-0.745 


-0.896 


SW-HSOW 


-0.685 


-0.474 


Ref. § 


-0.690 





FIG. 1. Transverse excitations u)k/J for 
the Kugel-Khomskii model at orbital degeneracy [Ez — 0) 
within RPA for the AFzz (top) and AFxx (bottom) phases 
in the fee (AFzz) Brillouin zone. Strong coupling between 
the (SW and SOW) modes results in a soft mode along the 
r - X(Y) direction. 



FIG. 2. The same as in Fig. |l|, but within the simplified 
SW-I-SOW scheme; the SOW dispersion is 1.5J. 



FIG. 3. Order parameters (Si) for AFzz (left) and AFxx 
(right) phases as functions of Ez/J using: full RPA (full 
lines), and SW-I-SWO scheme (dashed lines). The horizon- 
tal lines show the limits found a,t E^/ J —> -co (dashed line), 
and at E^/J oo (2D HAF, dashed-dotted line). 



FIG. 4. Ground state energies E of the AFzz (left) and 
AFxx (right) phases as functions of Ez/J, obtained using 
RPA (full hues), and SW-^SWO scheme (dashed lines). 
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